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Abstract 


A  technique  that  simultaneously  maximizes  the  gain  of  an  antenna  array  and 
places  a  number  of  independent  nulls  in  the  radiation  pattern  has  already  been 
presented.  The  elements  used  were  idealized  point  sources.-  In  the  present  paper, 
the  previous  analysis  is  extended  to  cover  arrays  of  thin  wire  elements  in  which 
interelement  mutual  coupling  and  scanning  effects  are  included.  Sample  calcula¬ 
tions  and  computed  radiation  patterns  demonstrate  that  an  appreciable  amount  of 
pattern  and  sidelobe  control  can  be  obtained  with  only  a  small  sacrifice  in  gain. 
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Null  Steering  arid  Maximum  Gain  in 
Electronically  Scanned  Dipoie  Arrays 


1.  INTRODUCTION 

In  a  recent  paper  (Drane  and  Mcllvenna,  1970),  the  authors  presented  a 
technique  for  maximizing  the  gain  of  ar.  antenna  array,  while  at  the  same  time 
placing  one  or  a  number  of  independent  nulls  and/o"  Celebes  in  the  far-field 
radiation  pattern.  *  Some  problem  f.reas  requiring  such  a  technique  include  the 
elimination  of  jammer  interference  in  radio-radar  commimioahnns  links,  the  re¬ 
duction  of  ground  reflections  in  sited  antennas  and  thr.  minimization  of  inter¬ 
antenna  interference  effects  in  muUi*an<enna  environments.  The  sample  computa¬ 
tions  in  that  earlier  paper  were  restricted  to  arrays  of  idealized  antenna  elements, 
that  is,  isotropic  radiators.  It  was  pointed  out  there,  however,  that  the  "maximum 
gain  -  null  placing  .jehnique  would  be  epp.icable  even  with  arrays  of  e,  omento 
(dipoles,  for  example)  when  in  ter  element  mutual  coupling  effects  must,  of  neces¬ 
sity,  be  considered.  Since  that  time,  a  particular  method  of  design  /or  arrays  of 
wire  elements  (Strait  and  Hirasawa,  1968  a.id  1968a)  which  elegantly  and  practically 
accounts  for  mutual  coupling  effects  has  been  combined  with  the  afore mentioned 
maximum  gain  -  constraint  technique.  These  design  methods  are  cib.eusr.ed  below 
together  with  some  examples  of  practical  interest. 

(Received  for  publication  31  January  1972) 

♦Several  papers  on  this  and  closely  related  subjects  have  appeared  (Nemit,  1969; 
Sandrin  and  Glatt,  1970;  Riegler  and  Compton,  1970;  Pierce,  1970;  Hesse)  and 
Sureau,  1971;  Adams  and  Strait,  1970;  Sanzgiri,  etal,  1971). 
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2.  DESIGN  TECHNIQUES  FOR  ARRAYS  OF  WIRE  ELEMENTS 

The  method  summarized  below  is  due  to  Strait  and  Harrington  and  has  been 
reported  in  detail  by  Strait  and  Adams  (1970),  Chao  and  Strait  (1971),  Strait  and 
Hirasawa  (1970a),  Hirasawa  and  Strait  (1971a,  1971b),  and  Strait  and  Chao  (1971). 
With  it,  the  designer  can  determine  such  properties  as  the  actual  current  distribu¬ 
tions  along  the  antenna  elements,  input  and  driving-point  impedances  and  the  radi¬ 
ation  and  scattered  field  patterns  of  arrays  of  parallel,  staggered  or  even  arbitrarily 
bent  wires.  The  elements  may  be  arbitrarily  spaced,  loaded  or  unloaded,  lossy  or 
lossless,  center-fed  or  arbitrarily  excited  at  points  along  their  length,  and  they 
may  be  placed  in  any  geometric  arrangement.  Even  radiating  ilements  with  wire 
junctions  such  as  crossed  dipoles  (Chao  and  Strait,  197 1)  can  be  handled.  Most 
importantly,  in  every  case  mutual  coupling  effects  are  included.  Detailed  handbook- 
typi  -jports  with  computer  programs  for  a  wide  variety  of  array  design  problems 
are  available  by  Strait  and  Hirasawa  (1969b,  1970b),  Chao  and  Strait  (1970),  and 
Mautz  and  Harrington  (1971).  Only  the  highlights  of  the  technique  are  presented 
here;  they  lead  to  an  expression  for  array  gain  with  mutual  coupling  effects  included. 

Starting  from  the  integral  equation  for  a  thin  wire,  one  considers  each  antenna 
in  the  array  to  be  subdivided  into  a  number  of  segments,  K  (ten  per  wavelength  has 
proven  to  be  sufficient  for  far-field  considerations),  V  ose  end-points  are  treated 
as  the  terminals  of  a  multi-terminal  pair  network.  The  goal  is  to  relate  every 
terminal  pair  to  every  other  terminal  pair  through  a  mutual  impedance  matrix. 

The  two  boundary  conditions,  that  the  tangential  component  of  the  electric  field 
vector  E  must  be  zero  on  the  surface  of  each  conducting  element  and  that  the  axial 
current  must  go  to  zero  at  the  ends  of  the  £  itenna,  are  included  in  the  analysis. 

Introducing  column  vectors  I  and  V,  each  with  K  components  that  are  respec¬ 
tively  the  segment  currents  and  voltages, 

T 

I  -  {  I^,  Ig,  •  •  • ,  ) 


and 


VT  =  {  Vj.  V2 . vK} . 

where  superscript  T  denotes  the  vector  transpose  operation,  allows  one  to  repre¬ 
sent  the  array  problem  in  terms  of  circuit  theory  format,  that  is, 

V  *  Z  I  , 
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where  Z,  a  (K  x  K)  -  element,  symmetric  complex  matrix,  is  the  mutual  impedance 
matrix  relating  the  voltage  and  current  for  any  segment  to  the  voltage  and  current 
for  every  other  segment. 

Methods  for  calculating  the  elements  in  Z  for  any  array  geometry  are  available 
(Strait  and  Hirasawa,  1969a).  It  has  been  shown  that 


Zk,nBjw'ioAinAik0(n*k) 


[tf/(n  ,k  )  -  0(n',k‘)  -  ^(n  ,k  )  +  0(n  ,k  )  ], 


where 


*<n'k>  =  4557 


irj 

n  J&t 


exp  [-jPlW2')) 

- Mi -  ri-i 


Rkn<z’> 


cz'. 


tti  th 

and  A*n  is  the  length  of  the  nm  segment,  Rkn(z')  is  the  distance  between  the  k 

and  n^‘  segments  and  the  superscripts  +  and  -  (in  n+  and  n"),  denote  the  starting 
point  and  termination,  respectively,  of  the  n**1  segment.  The  quantity  u>  is  the 
operating  angular  frequency,  /3  is  the  corresponding  wave  number,  and  e0  and  hq 
are  the  universal  electric  and  magnetic  constants.  The  only  other  approxima¬ 
tions  made  in  this  analysis  are  that  the  current  and  charge  are  assumed  to  be  con¬ 
stant  over  the  length  of  each  segment.  Calculations  and  comparisons  with  other 
mutual  coupling  methods  have  shown  that  such  assumptions  are  not  overly  restric¬ 
tive  (Strait  and  Hirasawa,  1969a).  Resistive  or  reactive  loading  and  element 
losses  can  be  included  simply  by  adding  a  load  matrix  Z{  to  the  Z  matrix  defined 
above  (Hirasawa  and  Strait,  1971a). 

Given  an  array  geometry  and  a  set  of  driving  voltages,  the  segment  currents, 
with  all  mutual  effects  included,  are  found  quite  simply  from 

I  =  Z*1  V  =  YV, 


where  Y  is  the  (KxK) admittance  matrix.  With  these  segment  currents,  all  the 
parameters  of  interest  to  the  array  designer  can  be  obtained.  For  example,  the 
far-zone  vector  potential  for  z-directed  current  elements  has  only  a  z-component 
given  by  (Strait  and  Hirasawa,  1967a) 
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and  the  far-field  amplitude  pattern  is 

E  (0,  4>)  =  -jw  »0‘‘-x  <®.  v)  sin  0. 

In  these  expressions  xQ,  yn  and  zn  are  the  Cartesian  coordinates  of  the  center  of 
the  segment  while  r,  0  and  $  denote  the  polar  coordinates  of  the  far-field 
observation  point. 

A  non-restrictive  assumption  that  simplifies  computations  is  to  take  all  seg-- 
ment  lengths  to  be  equal.  Defining  an  angle -dependent  row  vector  F(0,  <f>)  with  K 
elements,  that  is, 

Fr(0,  4)  =  exp  [  (xn  sin  0  coo  v  +  yn  sin  0  sin  <l>  +  zn  cos  0)  ] ,  1  <  n  <  K, 

and  replacing  I  by  YV  leads  to  a  matrix  form  for  Az,  viz. , 

A,  (0.  +  )  ■  Kn  [F(0,*)YV]t  (2) 

z  o 

where 


k  *  ‘-gEi:^r).  hi, 

o  4  »r 


where  A  is  the  length  common  to  all  segments.  One  need  only  use  the  relations 
above  ano  some  standard  definitions  to  develop  an  expression  for  antenna  gain 
which  includes  the  effects  of  mutual  coupling  (Strait  and  Hiraeawa,  1963a). 

Recall  that  power  gain  for  ar.  anter.ru  is  defined  to  be 


_  ..  power  radiated  in  a  particular  direction 
total'power  input  to  array 


(3) 


The  power  radiated  in  a  particular  direction  (6q,  *>0)  is  elated  to  the  far-field 
amplitude  pattern  by 


P  «V  V « 


r2/E  (0Q.  *Q>/2 


where  is  the  impedance  of  free  space.-  in  terms  of  the  matrix  expressions 
above. 
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/  uuAf  Y  _+  _  „ 

P<eo*^o)  “  n  \ — S — )  lV <F'Y)+  (FY)V)  sin- i  6, 

where  superscript  +  denotes  the  combined  operations  of  complex  conjugation  and 
transposition.  The  denominator  term  in  Eq.  (3)  can  be  handled  by  multi-port 
circuit  theory  techniques  and  is  expressible  in  terms  of  the  D  driving -point  voltages 
vd  and  currents  id  an  'Strait  and  Hirasawa,  1969a) 

D 

PT  =  E  Real  (Vd  ld*}  * 
d  =  1 

Let  and  VD  represent  D-element  vectors  of  the  driving-point  currents  and 
voltages.  In  vector  notation 


PT  *  Real  (VD  ID*). 

Let  Yd  represent  the  (D  x  D)  complex  driving-point  admittance  matrix  formed 
from  Y  by  selecting  only  those  elements  that  correspond  to  driving  points  (Strait 
and  Hirasawa,  1969a).  Because  of  the  symmetry  in  YD, 

pT  =  vDf  <Real  (yd,,vd 


and  Eq.  (3)  becomes 


G  =  K 


V  (FY)'  (F  V)  V 
1  VD+(ReaUYD))VD 


where 


„  (un  At) 

Ki  =  sin  6  -fer 


(4) 


—  +  — 

the  numerator  matrix  (F^)  (FY)  is  a  (K  X  K)-element  matrix  and  the  c'erominator 
matrix.  Real  (Yp),  is  a  (D  x  D)-element  matrix. 

In  the  authors'  earlier  paper  (Drane  'd  Mcllvenna,  1970),  the  parameter 
optimized  was  directive  gun  in  contrast  to  ihe  po^er  gain  discussed  above.  The 
difference  between  the  two  gains  lies  only  in  the  denominator  term  of  Eq.  (3). 
Directive  gain  utilizes  the  total  power  radiated,  that  is. 
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2»  v 


•Rad.  =  /  /P(0' 


<i>)  sin  8  d  8  d 


o  o 


instead  of  the  total  power  input  to  the  ar’-ay.  When  the  radiating  elements  are  loss¬ 
less  (the  usual  assumption),  the  two  gain  formulations  are  identical.  In  any  case, 
however,  directive  gain  is  an  upper  bound  for  the  power  gain  and  is  therefore  often 
used  as  a  comparative  indicator  of  array  performance.  It  can  be  calculated  from 
the  expressions  for  P(0, 4>)  given  earlier.  (See  Appendix  A.  for  details. ) 

Note  that,  as  yet,  no  assumptions  about  trie  number  of  feed  points  have  been 
made.  If  for  example,  every  segment  is  fed,  then  D  =  K,  Vq  =  V,  and  both  the 
numerator  and  denominator  matrices  in  Eq.  (4)  are  (K  x  K)  -element  matrices. 
Usually,  however,  D  <  K,  and,  in  fact,  a  very  common  array  configuration  has 
each  antenna  element  center-fed,  in  which  case,  D  =  N.  This  center-fed  assump¬ 
tion  also  provides  some  simplification  in  the  numerator  of  Eq.  (4).  For  perfectly 
conducting  center-fed  elements,  the  only  non-zero  entries  in  V  are  at  the  N  driving 
points  themselves.  One  may  then  delete  all  the  elements  in  the  numerator  matrix 
(FY)+  (FY)  which  do  not  correspond  to  these  feed  pcints,  denoting  the  resulting 
(N  x  N)  deleted  matrix  as:. 

A  =  (FY)d  (FY)d  . 


For  example,  in  a  two-element  array  with  five  segments  per  element  (K  =  10),  the 
center-feed  points  correspond  to  segments  3  and  C,  and  one  need  retain  only  four 
of  the  original  100,  that  is  (K  X  K)  elements  of  tne  numerator  matrix  (FY)  (FY), 
namely,  those  with  corresponding  pairs  of  indices  (3, 3),  (3,8),  (8,  3)  and  (8, 8). 
Thus,  the  center-fed  assumption  reduces  the  dimensions  of  the  numerator  and 
denominator  matrices  of  Eq.  (4)  from  (10  X  10),  that  is  (K  X  K)  to  (2  X  2),  that  is, 
(N  x  N).  The  final  form  for  the  power  gain  of  an  N-element  array  of  center-fed 
thin  wires,  with  mutual  effects  accounted  for,  can  then  be  wiitten  as 


G  =  K 


_  VD  A> 

1  VD  Real(YD)VD 


K, 


VD  *  Vp 
VD+SVD 


(5) 


The  numerator  and  denominator  matrices  in  Eq.  (5)  are  both  Hermitian,  A  is 
a  one-term  dyad,  B  is  positive  definite,  hence  the  maximum  gain  can  be  v/ritten 
down  immediately  as  (Strut  and  Hirasawa,  1969a) 


ginax 


-  K, 


(Real  (YD))_1(FY)D  , 


(6a) 


•U4  gxt 
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and  the  voltages  that  produce  it  are  given  (to  within  a  constant)  ay 

ytnax  r  (Real  (yd))_1(FY)d  .  (6b) 

Equations  (5)  and  (6)  serve  as  the  starting  points  for  the  maximum  gain-constraint 
technique . 


3.  THE  CONSTRAINT  METHOD 

The  technique  for  placing  nulls  in  the  far-field  radiation  pattern,  while  simul¬ 
taneously  maximizing  gain,  follows  the  steps  outlined  in  the  authors'  earlier  paper 
(Drane  and  Mclh  cnna,  1970).  Placing  nulls  in  the  directions  {  Z^,  <t>i  }, 
i  *  1,  2,  ....  M,  with  M  <  (N  -  1),  requires  that  E  (0, 4>)  (or  eoaivalently  Az(9,  <f>) ) 
be  set  to  zero  in  the  M  desired  directions.  For  the  center-fed  wires  discussed 
earlier,  we  therefore  generate  a  set  of  M  linear,  homogeneous  constraint  equations 
of  the  form 

(Fld^jlYlpVjj-O,  i  =  1,2,  ....  M.; 

Each  of  these  equations  consists  of  the  inner  product  of  a  D-element  angle- 
dependent  row  vector 

(F(0.,^)Y)D  . 

henceforth  referred  to  as  the  constraint  vector,  and  the  column  vector  of  driving 
point  voltages,  VD.  These  constraint  equations  are  not  directly  solved  for  the 
ci  mg  voltages.  Instead,  their  effect  is  introduced  into  the  g?.ir.  expression 
through  a  transformation  as  follows.,  Each  one  of  the  M  constraint  vectors  is 
considered  to  be  a  particular  row  in  a  constraint  matrix  which  thus  has  M  rows 
and  D  columns.  In  the  case  of  the  center-fed  wires  being  considered  here,  D  =  N. 
The  rectangular  (M  x  N)  constraint  matrix  is  made  into  a  square  (N  x  N)  matrix, 

C,  by  adding  (N  -  M)  rows.  These  additional  rows  can  be  any  collection  of  N- 
element  independent  vectors.  For  example,  the  simplest  such  collection  could  be 
(N  -  M)  rows  or  columns  of  the  (N  x  N)  identity  matrix..  The  (IJ  X  N)  matrix  C  is 
then  orthogonalized,  row  by  row,  using  any  appropriate  method  (Guillemin,  1949); 
and  the  resulting  (N  X  N)  matrix,  denoted  as  P,  contains  all  the  constraint  informa¬ 
tion  and  is  now  introduced  into  the  gain  expression  via  the  transformatic n 

VDr?+v  cr  V»P/D. 
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Thus,  from  Eq.  (5), 


G  =  K 


V+ACV 

1  +  — 

v  B  V 

c 


where 

Ac  =  P  A  P+  =  P  (FY)*  (F  Y)d  P+  =  acac+ 


and 


Bc  =  P  (Real  (Yp)  )  P+  . 


Following  the  arguments  outlined  earlier,  we  discard  the  first  M  rows  and  columns 
in  Ac  and  Bc>  leaving  the  abridged  (N  -  M)  x  (N  -  M)  element  matrices  Aftand  Ba< 
discard  the  first  M  entries  in  V  (denoting  the  resultant  abv  laged  vector  V .  ) 
and  finally  discard  the  first  M  entries  in  ac  (denoting  xhe  resultant  vector  aa).  Thus, 
the  final  quadratic  form  representation  for  array  gain,  with  all  the  effects  of 
mutual  coupling  and  pattern  constraints  included  is 


K 


V-  X. 


a  va 


1  V+B  V 
a  a  a 


the  maximum  constrained  gam  is 


pjnax 


K1  K  B."‘  »a> 


^  ^max 


and  the  corresponding  driving  voltages  are  found  from 


vr*  -  K  k1  v 


where  P^  is  an  (N  -  M)  x  N-element  matrix  formed  by  deletion  of  the  first  M  rows 
in  the  (N  y  N)-element  matrix  P., 

These  results  verify  the  claim  made  in  our  earlier  work  that,  even  when 
mutual  coupling  is  included,  the  form  of  the  gair  relation  and  the  properties  of  the 
numerator  and  denominator  matrices  are  unchanged.  * 


*An  alternate  constraint  technique  has  been  proposed  by  Adams  and  Strait  (1970). 


4.  SOME  EXAMPLES 


One  of  the  most  common  linear  array  configux  ^tions  utilizes  */2  elements 
uniformly  spaced  at  X/2  intervals.  This  geometry  was  used  for  the  sample  com¬ 
putations  in  this  paper.  There  were  12  thin,  straight,  wire  elements,  oriented  in 
the  z -direction  and  spaced  along  the  x-axis,  with  5  segments  per  element.  Strait 
and  Hirasawa  (1969b,  1970b),  Chao  arid  Strait  (1970),  and  Mautz  and  Harrington 
(1971)  provide  the  computer  programs  to  consider  other  geometries  and  element 
types. 

One  of  the  problems  ideally  suited  to  the  maximum  gain  -  constraint  technique 
is  the  elimination  of  interference,  jamming  or  unwanted  signals  through  control  of 
the  radiation  pattern.  Unconstrained  maximum  gain  radiation  patterns  are  charac¬ 
terized  by  the  presence  of  a  relatively  high  secondary  peak,  for  example  13  dB 
below  the  main  peak.  These  high  secondary  peaks  make  such  arrays  susceptible 
to  jamming.  One  can  use  constraints  in  several  different  ways  to  reduce  these  high 
peaks;  two  possibilities  are  discussed  in  the  examples  below. 

One  corrective  technique  requires  the  use  of  a  single  constraint  and  replaces  a 
secondary  peak  of  the  unconstrained  maximum  gain  pattern  with  a  null.  A  second 
approach  is  to  use  two  closely  spaced  nulls  to  reduce  the  radiation  level  over  an 
angular  sector  that  includes  the  secondary  peak.  The  first  approach  is  useful  when 
the  angular  location  of  the  jamming  source  is  exactly  known.,  The  second  is  more 
suitable  when  the  location  of  the  interfering  source  is  not  accurately  known  or  when 
one  desires  to  control  a  broader  region  in  the  pattern  perhaps  because  the  jamming 
source  is  of  some  s;<"iificant  angular  extent.  In  both  cases,  the  null  or  nulls 
should  be  maintained  fixed  in  space  even  while  the  main  beam  is  s*.  anned  all  the 
way  from  broadside  to  endfire.  Because  of  this  wide  range  in  scan  angles,  it  is 
especially  important  to  include  mutual  coupling  effects..  Figure  1  demonstrates 
the  two  techniques  with  the  unconstrained  broadside  pattern  of  a  12 -element,  X/2  - 
spaced  array  of  thin-wire  dipole  elements,  ty2  long.  All  patterns  are  normalized 
to  facilitate  comparison.  The  secondary  peak  of  the  unconstrained  maximum  gam 
pattern  at  0  =  '’G0  is  eliminated  by  the  first  technique  or  it  can  be  reduced  by  the 
second  technique.;  Not'*  that  the  pattern  structure  is  relatively  unchanged,  except 
in  the  vicinity  of  the  enforced  nulls.;  This  demonstrates  the  localized  pattern 
control  available  with  the  constraint  method.  The  3  dB  beamwidths  of  all  three 
patterns  are  almost  identical,  while  the  beamwidth  between  nulls  of  the  constrained 
patterns  increases  only  slightly.  For  the  two -constraint  case,  the  design  criteria 
were  as  follows.  The  main  beam  should  be  scannable  in  the  principal  H-plane,  but 
for  every  scan  angle  in  that  plane,  the  gain  in  the  direction  of  the  scan  angle  is 
maximum,  while  the  radiation  level  is  held  at  or  below  -30  dB  over  a  4°  angular 
sector  centered  at  <t>  =  76°.  This  minimal  radiation  level  sector  was  created  by 
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Figure  1.  Broadside  Beams 


using  two  nulls,  fixed  in  place  (at  6  =  74°  and  <)>  =  78°)  for  all  scan  angles.  The 
resultant  patterns,  for  this  two-constraint  case  plotted  at  10°  increments  in  the 
scan  angle,  are  shown  in  Figure  2.  Note  in  Figure  2(a)  the  well  defined  trough 
in  the  patterns,  which  is  at  least  4°  wide  at  every  scan  angle  and  in  which  the  highest 
radiation  level  is  -30  dB  for  the  broadside  pattern  and  at  least  -35  dB  for  the 
other  patterns.  Figure  2 <b>  shows  the  gradual  and  expected  broadening  of  the  main 
beam  and  the  buildup  of  the  endfire  lobes  as  scan  angle  increases,  until  at  (/>  =  180°, 
the  pattern  is  bifurcated.  This  splitting  of  the  main  beam  into  two  equal  lobes  is 
characteristic  of  the  X/2  spacing  chosen  for  this  example.. 

Complete  information  for  the  unconstrained,  one -constraint  and  two-constraint 
cases  is  presented  in  Taole  1.  Corresponding  radiation  pattern  plots  are  shown  in 
Figures  1,  and  3  through  11. 

Examining  Table  1(a),  we  see  that  as  mentioned  earlier,  constrained  gains  are 
less  than  or  equal  to  the  unconstrained  gains.  An  interesting  observation  is  the 
small  amounts  of  ge‘-'  that  must  be  surrendered  to  obtain  a  desired  degree  of 
pattern  control!.  In  the  one-constraint  case,  the  largest  loss  occurs  at  broadside 
and  represents  only  a  5  percent  reduction  in  gain.  In  the  two -constraint  case,  the 
largest  loss  is  about  4  percent  again  at  broadside,  and  in  both  the  one-  and  two- 
constraint  cases,  the  jain  loss  is  less  than  about  2  percent  over  the  range  of  other 
scan  angles.- 


Table  1  (Contd).  Comparison  of  Unconstrained  (U),  One  Constraint  (Cl)  and  Two  Constraint  (C2)  Array  Properties 


UNCONSTRAINED 

1  CONSTRAINT 

2  CONSTRAINTS 


+  (d«gr«ts! 


Figure  3.  Beam  Scanned  to  100 


-UNCONSTRAINED 
- 1  CONSTRAINT 
•2  CONSTRAINTS 


4  (dtQrets) 


Figure  4.  Beam  Scanned  to  110° 


4 ( dtgrtts ) 


^  ( degrees) 


Figure  10.  Beam  Scanned  to  170° 
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Figure  11.  Beam  Scanned  to  180° 


We  also  see  that  for  some  scan  angles,  that  is  90°  to  110°,  the  gain  in  the 
two-constraint  case  i3  higher  than  in  the  one -constraint  case.  This  is  simply 
explained.  Controlling  the  value  of  the  secondary  peak  to  -30  dB  as  done  in  the 
two  constraint  case  represents  less  of  a  constraint  on  the  pattern  behavior  at  these 
angles  than  does  forcing  that  same  peak  value  to  be  zero,  as  is  done  in  the  one- 
constraint  case.  It  thus  does  not  necessarily  follow  that  the  gain  decrease  is 
directly  proportional  to  the  number  of  null  constraints  used.  Only  when  the  con¬ 
straints  are  used  to  perform  the  same  function  on  the  pattern  will  the  gain  decrease 
be  proportional  to  the  numb*,'  of  constraints.  For  example,  replacing  two  peaks 
with  nulls  causes  more  gain  loss  than  replacing  one  peak  with  a  null,  and  replacing 
three  peaks  with  nulls  causes  even  more  of  a  gain  loss,  and  so  on.  This  type  ol 
dependence,  linking  gain  loss  to  number  of  enforced  nulls,  is  the  only  one  that  can 
be  identified.. 

Generally  speaking,  the  decrease  in  gain  value  caused  by  pattern  constraints 
is  proportional  to  the  value  of  the  unconstrained  pattern  in  the  immediate  vicinity 
of  the  enforced  nulls.  For  example,  as  seen  above  and  in  Figure  1,  at  broadside 
the  constraint  legion  includes  a  significant  secondary  peak  of  the  unconstrained 
patter’  and  the  gain  loss  is  largest;  when  the  unconstrained  pattern  value  in  the 
constraint  region  becomes  smaller,  the  gain  loss  is  less.-  Fox-'  example,  at  140° 
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and  160°,  the  unconstrained  pattern  already  ha3  a  null  at  $  =  76°,  implying  that  the 
one  constraint  case  is  in  actuality  no  constraint  at  all.  There  is  therefore  no  loss 
in  gain,  the  amplitude  and  phase  distributions  are  almost  identical  and  the  uncon¬ 
strained  and  one  constraint  patterns  are  quite  similar  (see  Figures  7  and  9). 

An  examination  of  the  voltage  amplitude  distributions  in  Table  1(b),  (each 
distribution  is  normalized  to  the  voltage  for  its  first  element)  reveals  that  these 
maximum  gain  voltages,  constrained  and  unconstrained,  are  all  non-uniform  and 
most  are  non-eymmetric  about  thfe  center  of  the  array.  The  only  two  symmetric 
distributions  (the  broadside  and  endfire  unconstrained  cases)  correspond  to  the  only 
two  radiation  patterns  that  are  symmetric  about  the  broadside  angle  <f>  =  90°.  Note 
that  in  some  cases  the  voltage  required  for  element  12  differs  significantly  from 
the  voltages  for  all  other  elements.  It  also  appears  that  for  certain  scan  angles, 
for  example  150°,  the  difference  between  the  unconstrained  and  constrained  ampli¬ 
tude  distributions  is  small.  As  was  true  for  the  earlier  comparison  of  the  corre¬ 
sponding  gain  values,  the  explanation  here  is  that  at  some  scan  angles  the  uncon¬ 
strained  pattern  values  are  either  very  small  in  the  constraint  region,  or  the 
direction  for  one  of  the  unconstrained  nulls  is  nearly  the  same  as  that  of  one  of  the 
constraints..  The  constraint  condition  therefore  requires  only  a  slight  perturbation 
of  the  unconstrained  voltage  distribution  in  these  cases.  Finally,  note  that  in  no 
case  is  the  dynamic  range  of  the  driving-point  voltages  excessive.  Compared  to 
the  ten-to-one  variation  usually  taken  as  the  practical  limit  of  realizability,  the 
range  for  both  constrained  and  unconstrained  optimum  distributions  is  less  than 
about  two-to-one,  and  the  distributions  should  offer  no  construction  problems.  The 
phase  distributions  for  these  maximum  gain  arrays  can  be  considered  to  consist  of 
the  superposition  of  several  components.  The  first  is  the  uniform  progressive 
phase  taper  that  steers  the  beam  to  the  desired  scan  angle  and  is  common  to  all 
electronically  scanned  arrays,  const-  ained  or  unconstrained.  The  second  taper  is 
that  required  to  maximize  the  gain.  And,  in  the  constrained  cases,  an  additional 
taper  is  required  to  place  the  pattern  nulls  as  prescribed.  Table  1(c)  shows  the 
unconstrained  and  constrained  normalized  phase  distributions,  after  the  beam 
steering  phase  distribution  (shown  for  reference  in  Table  1(d))  has  been  removed. 

As  in  the  case  of  the  voltage  distributions,  note  here  also  that  the  phase  distribu¬ 
tions  are  non-uniform  and  that  most  are  non-symmetric.  Finally,  keep  m  mind 
that  the  examples  discussed  in  this  paper,  and  indeed  in  most  expected  applications, 
even  though  the  designer  directly  controls  the  sidelobe  level  in  only  some  local 
region  of  the  radiation  pattern,  additional  control,  as  a  direct  result  of  the  applica¬ 
tion  of  these  constraints  will  not  become  necessary  over  the  remainder  o"  the 
sidelobe  region  (Sanzgiri  and  Butler,  1971).  This  is  due  to  the  general  phenomenon 
that  reasonably  low  sidelobes  are  an  intrinsic  concomitant  to  maximizat'on  of  gain. 
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5.  CONCLUSIONS 

We  conclude  that  even  with  mutual  coupling  accounted  for,  the  maximum  gain  - 
constraint  technique  provides  realistic  and  practical  solutions  to  array  design 
problems  requiring  localized  pattern  control.  The  technique  is  applicable  to  prac¬ 
tically  all  arrays  of  wire  elements,  regardless  of  their  geometric  arrangement  a,  d 
their  electrical  loading  characteristics.  The  required  relations  and  equations  can 
be  cast  in  matrix  form  and  are  ideally  suited  to  computer  manipulation,  making 
them  useful  in  adaptive  or  real-time  situations  requiring  rapid  pattern  reconfigura¬ 
tion. 
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Appendix  A 

Derivation  of  a  Rotation  for  the  Directive  Gain  of  an  Array  of  Electrically 
Thin,  Straight  Wires  With  Mutual  Coupling  Effects  Included 


It  was  pointed  out  in  the  text  that  the  expressions  for  directive  gain  and  power 
gain  differ  only  in  the  form  of  the  denominator  term.  Power  gain  used  the  total 
power  input  to  the  array  while  directive  gain  uses  the  total  power  radiated.  These 
are  the  same  when  the  system  is  lossless.  In  all  other  cases,  directive  gain  serves 
as  an  upper  bound  on  the  power  gain  and  is  often  used  as  a  performance  criterion  in 
array  design.  This  latter  expression  is  given  in  general  by 


PRad. 
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2»  ir 
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where  P (0,<l>),  the  power  pattern  of  the  array,  can  be  expressed  as 


P (M) 
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E<9,  *) 


=  K 


A 2  sin2  6 


and  K  =  (rupo)2/no.-  This  integral  has  been  evaluated  by  Tai  (1964)  and  Lo  et  al 
(1966)  for  isotropic  sources  and  short  dipoles.  For  the  case  of  electrically  thin, 
z-directed  wires,  spaced  along  the  x-axis,  Az  (9,  <t>)  is  given  by 
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A z(9,+)  •  exP^jfir)  ^  AlnInexp  [  jp(xR  sin  flcos  ♦+ zn  cos  flj. 
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Taking  all  segments  to  be  of  equal  length  leads  to 
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and 


D  =  0  (x  -  x  ). 
nm  n  m 


E _ =  0  (z  -z  )  . 

nm  n  m 


The  b  terms  depend  only  on  the  array  geometry  and  in  no  way  are  modified  by 
the  inclusion  or  exclusion  of  mutual  coupling  effects.  But,  the  In  quantities  do 
indeed  depend  for  their  values  on  mutual  coupling  considerations.  In  matrix  nota¬ 
tion,  Eq.  (Al)  becomes 


p 

Rad 


V+  (Y+  0  Y)  V. 


(A2) 


The  integration  on  <t>  in  the  expression  for  bnm  is  readily  performed  (Jahnke 
and  Emde,  Pg.  149),  that  is. 
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This  integral  may  be  evaluated  as  follows.  Note  that  about  0  =  jt/2,  sin  0  and 
Jq  (Dnm  sin  0)  are  even  functions  of  0,  while  sin  (Enm  cos  0)  is  odd.  Hence,  we 
can  discard  the  imaginary  part  of  bnm  and,  using  the  evenness  of  the  real  part,  we 
have 
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(2i):  ■  r(l+2i)  =  r  [  2  (i+  1/2)] 


28 


and  that 
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and  b  =  2/3. 
nn 

Note  that 


:-»1+1/2'ii+,/2®m„> 

:-.)1+3/2jtt3/2 1.1, 

hence  b  =  b  and  8  is  a  real,  symmetric  matrix.  Both  m  and  n  range  from 
mn  nm  * 

1  to  K,  where  K  is  the  total  number  of  segments  in  an  N- element  array.  K  can  be 
rather  large;  for  example,  in  the  arrays  considered  in  this  paper,  12 -elements 
with  5  segments  per  element,  K  *  60.  Fortunately,  however,  one  needs  to  calcu¬ 
late  all  1800  of  the  bmn  elements  only  for  the  most  general  array  geometries,  that 
is,  wires  of  arbitrary  length  and  shape,  arbitrarily  spaced  and  located.  In  most 
arrays  of  practical  interest,  significantly  less  calculation  is  required.  For  equally 
spaced  arrays  of  N  straight  wires,  with  S  segments  per  wire,  only  NS  or  K  of  the 
1800  elements  in  0  are  distinct.  To  see  this,  refer  to  Figure  A1  taking  N  =  12 
elements  and  S  =  5  segments  per  element.  Each  of  the  twelve  antenna  elements 
gives  rise  to  25  bmn  terms.  These  (5  X  5)  -  element  square  submatrices  are 
situated  symmetrically  about  the  main  diagonal  in  0;  hence  the  srbmatrix  diagonal 
terms  are  known  to  be  2/3.  Because  of  the  symmetry  in  0,  only  10  of  the  remaining 
20  terms  in  any  main  diagonal  submatrix  can  be  distinct.  One  such  submatrix,  that 
associated  with  the  first  element,  is  shown  below. 


D  =  -  D 
mn  nm 


E  *  -  E 
mn  rii  • 


JUl/2®mn|si 
Ji+  3/2  ^mn*  = 


Each  one  of  the  10  required  b  terms  depends  on  both  a  D  and  an  E  term. 

nm  nm  nrn 

If  the  wire  is  straight  and  parallel  to  the  z-axis,  all  5  segments  on  that  wire  have 
the  same  x-coordinate  and  Dnm  «  0  for  all  10  of  the  desired  elements  >n  jj If,  ii 
addition,  all  segments  are  chosen  to  be  equal,  then 


Figure  Al.  12-Element  Array,  5 -Segments  per  Element 
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Hence, 

b12  *  b23  *  b34  *  b45 
b13  *  b24  “  b35 
b14  *  b25  * 

and  only  5  elements  are  needed  to  completely  specify  Its  final  form  is 


2/3 

b12 

b13 

b14 

b15 

b12 

2/3 

b12 

b13 

b14 

bl3 

b12 

2/3 

b12 

b13 

b14 

b13 

b12 

2/3 

b12 

b15 

b14 

b13 

b12 

2/3 

If  all  the  elements  in  the  array  are  identical,  then  all  of  the  12  main  diagonal  sub- 

matrJces  ar_-  identical  to  and  300  elements  in  P  are  determined  with  only  four 

b  „  calculations, 
mi 

The  off-diagcnal  submatrices  in  p  depend  on  the  array  spacing.  Each  pair  of 

elements  in  the  array  gives  rise  to  a  25  term  oubmatrix  and  for  a  K-element  array, 

there  are  66  of  these.  In  any  one  of  these  submatrices,  if  the  uvo  antenna  elements 

in  question  are  straight  and  parallel,  all  D  terms  are  the  same  and  have  the 

value  Pd  where  d  is  the  interelement  separation.  As  in  ^  there  are  only  5  distinct 

E  „  values,  hence  once  again  only  5  calculations  are  needed  to  determine  all  25 
nm  _ 

terms  in  any  of  the  submatrices.  The  submatrix,  representative  of  the  first 
two  array  elements,  is  shown  below.  It  has  the  same  type  of  symmetry  and  row 
regularity  as  py  sometimes  called  a  Toeplitz  property. 


b16 

b17 

00 

pH 

-O 

b19 

bl,  10 

b17 

b16 

b17 

b18 

b19 

b18 

b17 

b16 

bl7 

b18 

b19 

b18 

b17 

b16 

b17 

bl,  10 

b19 

b18 

b17 

b16 
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Similar  arguments  show  that  the  submatrices  representative  of  the  remaining 
array  elements  and  element  one,  that  is  0^,  04,  ,  0j2*  the  same  form, 

symmetry  and  regiu  rity  as  02  and  that  there  are  only  5  distinct  terms  in  each. 

When  the  array  is  equally  spaced,  each  of  the  remaining  48  submatrices  in 
0  is  identical  to  one  of  the  02.  0g.  ...  p..,  submatrices  discussed  above.  Thus, 
for  an  equally  spaced  array  of  straight  parallel  wires  with  equal  segments,  p  con¬ 
sists  of  only  12  distinct  submatrices,  and  in  these,  only  60  terms  are  distinct;  it 
has  the  same  Toeplitz  property  as  its  component  submatrices. 

^3  h  h  h  h  •  •  •  ^12 

^2  ^1  ^2  ^3  ^6  h  '  ’ '  ^11 

F3  h  h  h  ^3  Ms  P6  •  *  *  ho 

P  * 


^12^11^13^8^7  ^6  h  ’ '  ‘ 

The  cyclic  property  oi  f  and  its  submatrices  is  also  present  in  the  generalized 

impedance  matrix  t.  for  this  type  of  array  (Strait  and  Hirasawa,  1969a),  since 

elements  in  both  p  and  Z  depend  only  on  intersegment  distances. 

The  60  distinct  elements  in  0  are  the  first  row  blm>  1  £  m  S  60,  or  the  first 

column  terms.  These  60  terms  are  associated  with  four  combinations  of  D  __  and 

nm 

E  : 
nm 


(a) 

D 

=  E  =0 

nm 

nm 

(b) 

D 

0,  E  =0 

nm 

'  nm 

(c) 

D 

=  0,  E  i  0 

nm 

(d) 

D 

nm 

4  0,  E  f  0  . 

1  nm  ' 

These  can  be  used  to  simplify  the  bnm  expression  in  Eq.  (A4).  Case  (a)  has  been 

discussed  earlier  and  constitutes  the  main  diagonal  terms  in  0  where  bnn  =2/3. 

Case  (b)  is  the  situation  in  11  of  the  60  required  b  elements.  They  are  b. 

nm  J  la 

where  a  =  6,  11,  16, .  . ,  56. 


32 


Referring  to  the  general  bQm  expression  in  Eq.  (A4),  note  that  when  E^^  *  0,  all 
but  the  first  (i  -  0)  term  disappears  and 


J i/2  <Pto> 


J3/2  (Dto) 

V5- 


for  these  11  terms.  But. 


2D,  sin  D, 

J./2  'J  -nr2-  sf 


and 


J3/2  03 


In'  y  » 


sin  D 


la 


cos  D 


la 


°v 


D 


la 


Therefore, 


*  <Dla>2 


sin  D,  , 

-DT^  <-D  +  cos 

ut 


(A5) 


Note  that  Case  (b),  with  D _ +  0  and  «  0  corresponds  to  elements  of  zero 

length,  that  is  isotropic  or  short  dipoles,  and  has  been  treated  by  Tai  (1964)  and 
Lo  et  al  (1966).  Equation  (A5)  agrees  with  their  results.  Case  (c),  Dnm  *  0, 

E  _  4*  0  arises  in  calculating  4  of  the  required  60  b  „  elements.  They  are 
b^,  0*2, 3,4, 5.  Looking  at  a  general  series  representation  for  Jv  (x),  that  is. 


Jv  (x)  = 


y  (-Dr  xv+2r 

rl  2v+2r  r(v  +  r  +  1) 


we  see  that 


Jv  (x)  x 

— ~  -  — :  asx  °- 

xv  2  r(v+i) 

Hence,  if  in  Eq.  (A4)  D  =  0,  we  have 
nm 


blp 


•I 

1=0 


(-1)1  (E16)21 
(2i+l)!  (1  +  3/2) 


(A  6) 
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where  we  have  used  the  fact  that 


*1/2 

22l+1r(i+  1/2)  r(i-H) 


1 

2f'(2i  + 1)  - 


The  remaining  44  of  the  60  required  bnm  terms  must  be  calculated  from  Eq.  \A4). 

It  is  probably  easier  to  simply  use  Eq.  (A4)  for  all  60  b^  terms  when  programming 
a  computer  for  the  calculations.  In  this  case,  fifteen  terms  of  the  series,  Eq  (A4), 
will  provide  sufficient  accuracy,  that  is  on  the  order  of  10  ,  as  long  as  the  ir.ter- 

element  spacing  is  no  smaller  than  X/3.  Calculations  based  on  these  results  agreed 
with  those  based  on  the  power  gain  formulation  for  the  12 -element,  lossless  arrays 
considered  in  this  paper. 

Equations  (A2)  and  (A4),  togethei  with  the  symmetry  arguments  and  simplified 

b  forms,  provide  the  basis  for  calculating  the  directive  gain  of  an  array  of 
nm 

equally  spaced  straight  wires  with  mutual  coupling  effects  included.  The  finat 
formula  for  this  gain  is  given  by:. 
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